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ABSTRACT 


A  sensitivity  theorem  by  R,  M.  Cohn  states  that  for  linear  resistive 
circuits  the  ratio  of  a  fractional  change  in  the  d.c.  input  resistance 
to  a  fractional  change  in  an  internal  resistance  is  equal  to  the  square 
of  the  ratio  of  the  current  through  that  internal  resistance  to  the  d.c. 
input  current.   The  theorem  is  extended  to  show  the  sensitivity  of 
input  impedance  to  changes  in  internal  impedances  for  an  arbitrary  network 
at  all  frequencies.   Equations  are  developed  which  show  the  relation 
between  sensitivities  and  instantaneous  power  in  the  frequency  domain. 
An  extension  to  the  time  domain  makes  digital  computer  solutions 
possible . 
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I.   INTRODUCTION 

In  1950,  R.  M.  Cohn  [Ref.  1]  presented  a  linear  resistive  network 
sensitivity  theorem  which  was  new  and  significant  in  its  approach  to 
sensitivity.   Although  the  theorem  is  restrictive  in  its  original  form, 
it  is,  in  fact,  a  particular  case  of  a  more  general  result  involving 
circuit  sensitivities  as  shall  be  discussed  here. 

In  Chapter  2,  a  proof  of  Cohn's  Theorem  is  presented.   In  Chapter  3, 
the  theorem  is  extended  to  an  arbitrary  input  impedance.   In  Chapter  4, 
the  results  of  the  previous  chapter,  which  were  developed  in  the  frequency 
domain,  are  interpreted  in  terms  of  time  domain  responses.   A  general 
relationship  between  sensitivity  and  power  relationships  is  developed. 
The  results  are  then  illustrated  by  means  of  an  example  for  which  a 
computer  program  has  been  written  to  obtain   the  appropriate  waveforms. 


II.   PROOF  OF  COHN'S  THEOREM 

Consider  a  two  terminal  resistive  network  consisting  of  branches 

b-±   ,b2  , b  ,  having  indepenclendent  variable  resistors  ,  rj  ,r2  ,  .  .  .  .  ,r  . 

The  input  resistance  R  (see  Figure  1)  is  then  a  rational  function  of  the 
r's.   Consider  also,  that  a  unit  current  is  applied  to  the  input  terminals. 
It  can  then  be  shown  that  R  has  the  following  property: 

.8       P 

i  r  r 

BR     .  2       ,   oR  /R     mm  m                  , ,  N 

—  =  x      and  -r j —  =  — - —  =  — —                   (1) 

3r     m        or  /r      R  Pn 

m               mm  K 

where  i   is  the  current  flowing  in  b  .   Equation  (1),  which  is  Cohn's 

m  m 

Theorem  [Ref.  l~]  >    states  that  the  ratio  of  the  fractional  change  in  the 
input  resistance  to  the  fractional  change  in  an  internal  resistance  is 
proportional  to  the  ratio  of  the  power  dissipated  in  the  internal  resistance 
to  the  input  power.   The  proof  that  follows  is  based  on  that  of  Cohn. 

Let  i.  denote  the  current  flowing  in,  and  v.  the  voltage  across,  branch 
b..   Let  the  network  co.isist  of  q  nodes  designated  by  nx  ,n2  ,  .  .  .  ,n   with  the 
current  source  connected  from  nj  to  n  .   If  the  loops  are  designated  by 


h  »l2  »  •  •  •  >lt  an 


d  the  direction  of  current  flow,  in  each  branch  and  in  each 


loop  assigned  arbitrarily,  then  Kirchoff's  equations  may  be  stated  as: 
t 

)  b.  .v.  =  0        i  =  1,2, ,t  (2) 

^  ij  J  • 


)  a.  ,i.  =  g.       i  =  1,2, ,q  (3) 

-J  ij  J     J 

J-l 

b. .  =  1  if  the  branch  current  and  loop  current  have  the  same  direction, 

-   -1  if  the  branch  current  and  loop  current  have  opposite  directions, 
=  0  if  the  branch  is  not  in  the  loop. 
a.  .  =  1  if  the  branch  current  is  directed  away  from  the  node, 
=  -1  if  the  branch  current  is  directed  toward  the  node, 


=  0  if  the  branch  is  not  connected  to  the  node. 

g  -   1  for  j  =  1 

-  -1  for  j  ■  q 

=  0  elsewhere. 

Define  e.  =  i.r   to  be  the  voltage  across  b .  in  a  positive  direction 
J     J  J  J 

and  let  the  voltage  of  node  n.  with  respect  to  node  n,  be  E,  ,  =  f  e  with 

1  J      ij     J  j 

f.  =  1  if  i.  goes  from  n.  to  n  .  and  f,  =  -1  if  i   goes  from  n   to  n..   If 
J         J  i      J       J  J  j      i 


node  q  is  grounded,  E,   =  E,  =  R 

lq     1 


(4) 


Since  Kirchoff's  laws  imply  the  conservation  of  power,  the  power  input 
must  equal  the  sum  of  the  powers  consumed  by  the  individual  branches. 
This  may  be  written  as: 


n 

V  2 
E  I  =  )  i  .r  .  =  R 

1    ^  J  J 
J-l 


(5) 


If  equation  (3)  is  differentiated  with  respect  to  r  ,  it  is  found 

m 


that 


n 


Sr 
m 


,  1 


LLaij   jj 
J-l 


di  .   dg. 


=0    i  =  1,2, .. . ,q   (6) 


J-l 


ii  3r    3r 
m     m 


Let  the  nodes  of  b  .  be  denoted  by  n   and  n  ,  such  that  e   is  positive 

J  ja      jb       •     j 

going  from  n.   to  n.,  ,  then 
ja     Jb 


i.r.  =  e.i.  =  (E.,  -  E.  )i 
J  J     J  J    '  Jb    ja   j 


(7) 


n        n 

)    iZr.  =  }  (E.,  -  E.  )i. 
'-  J  J   Lj        jb    ja   j 


(8) 


J-l       3=1 


By  defining  a  function  h.  (an  arbitrary  function  of  the  r's)  such  that 


n 


;  a  h.  =  0 

-  ij  J 
J-l 


i  =  1,2, ,q 


(9) 


and  using  equation  (7)  above, 

n  n        n 

)  i.r.h.  =  /  e.h.  -  /(E   -  E.  )h.     .  (10) 

{->     J  J  J   *-.     j  j   L>     jb    ja   j 

j=l        j=l       j=l 

Since  E  .,  and  E.   are  voltages  measured  with  respect  to  node  q,  they 
jb       ja  r  n'     J 

may  be  written  as  E.'s.   Also,  a  network  is  not  necessarily  a  series 

connection  of  resistors.   Thus,  E  .,  and  E.   may  be  associated  with  more 

Jb       ja 

than  one  branch.   If  the  sum  on  the  right  hand  side  of  equation  (8)  is 
performed  and  all  terms  with  the  same  E.  are  combined,  it  is  seen  that: 

n  n  q   n 

)  i.r.h.  =Y  (E   -  E   )  h.  =Y  (  Ya.  .h  .  V         (11) 
L>     j  j  j   L       jb    ja   j  L»  \  U  lj  j/  x 

j=l        j=l  i=l  j=l 

=  0  by  the  definition  of  }    a. ,h . 

^  ij  J 

This  can  easily  be  shown  in  a  simple  example  (see  Figure  2). 

q 

I  (\b   "    Eka)hk   "    (E1    "   V"l   +   (E1   -    E3>"2    +   (E1   "   Vh5   + 
k=l 

(E2    -   E3)h4   +   (E3    -   Eq)h3 

=  E1(h1  +  h2   +  h5)    +  E2(h4   -   h5)    +  E3(-h2  +  h3   -   h4)   + 

E„(-h1    -   hj        .  (12) 

q         1  3 

From    the   definition    for   a. .    given   earlier: 

all   =  al2   =  a!5   =  a24   =  a33   =    l* 

a25  =  a34  =  332  =  aql  =   aq3  =  '1» 
and  all  others  equal  0. 

If  each  h   in  equation  (12)  is  multiplied  by  its  appropriate  a. . 
such  that  i  corresponds  to  the  appropriate  index  of  the  E. ,  it  will  be 
seen  that  equation  (11)  holds. 

Comparing  equation  (6)  and  the  defining  equation  for  h.,  (9),  it  is 

seen  that  di./dr   =  h  .   Thus  from  (11) 
J    m     j 


n  n      _  . 

Bi  . 


V  V 


J 


)  i.r  .h.  =  }    r  .i  .  — - J  =  0  .                    (13) 
(->     3    3    3        ^  ]  J  cr 

J-l        J-l 

If  equation  (5)  is  now  differentiated  with  respect  to  r  ,  the  result  is: 

m 

n          n    *  n      *  ■ 

Br  .  — ,     dx 


c>r        b    r  v  -2    1     V  -2_i  j.  o  V  •       c_j 


dr    Br   L  zL  jr jj   L     j  Br 

m     m    .  ,        .  ,     m      .  ,       m 


From  equation  (13) 

Bi 


/  l.i 

t-  J  J  orr 

j-l 


—J  =  0 
m 


(14) 


Since  the  resistors  of  the  circuit  are  independent, 

Br     |  0,  j/m 

3r     \  1 ,  j  =  m 
m 

and  equation  (14)  becomes 

^  -  i2     •  (14a) 

dr     m 

m 

This  is  Cohn's  theorem  for  a  resistive  network  with  a  unit  step  input 
current . 


III.   EXTENSION  OF  COIN'S  THEOREM  TO  AN  ARBITRARY  IMPEDANCE 

The  expression  developed  in  equation  (14a)  is  interesting,  but  is 

limited  to  purely  resistive  circuits.   Therefore,  an  extension  of  equation 

(14a)  to  an  arbitrary  impedance  is  desireable.   By  applying  a  theorem  of 
Tellegen  [Ref.  2],  it  can  be  shown  that: 

msi  ^[Vilf  (15) 


BZb(s)    Ll(s) 

The  method  used  follows  to  some  extent  the  work  of  Penfield  [Ref.  3]. 
Tellegen's  theorem  states,  "In  a  network  configuration,  imagine 
branch  currents  i  such  that  for  every  node  Si  =  0 ,  imagine  branch  voltages 
v  such  that  for  every  mesh  £v  =  0 ,  and  for  every  branch  let  the  positive 
direction  of  current  be  from  the  +  to  the  -  denoting  the  positive  polarity 
of  the  voltage.   Then  £iv  =  0,  where  the  summation  is  over  all  branches". 
Tellegen  then  goes  on  to  show  how  for  Figure  (1)  the  theorem  can  be 
written  as  £i  v   =  £  i.  v  where  i   and  v   are  the  terminal  voltages  and 

t  t  t     T.    D  D  t         t 

currents  and  i,  and  v   are  branch  currents  and  voltages.   It  is  significant 
b      b 

to  note  that  Tellegen's  Theorem  is  strictly  topological  and  is  related  in 
no  way  to  the  linearity  or  nonlinearity  of  the  circuit  components.   For  a 
constant,  linear  system,  Tellegen  shows  that  £l  V '  =  Sl/V'  .   The  unprimed 
quantities  are  the  first  state  values  and  the  primed  quantities  are  the 
second  state  values  of  two  circuits  with  the  same  topology.1 

Since  the  Laplace  operator  is  a  linear  operator,  it  can  be  applied 
and  Tellegen's  theorem  can  be  written  as: 


1  It  is  unfortunate  that  Penfield  chose  to  use  the  term  "state  values" 
since  it  can  easily  be  confused  with  the  "state  variables"  often  used  in 
circuit  analysis.   As  used  here,  "state  values"  refer  to  the  values  of  the 
voltages  and  currents  of  the  branches  of  two  different  circuits.   These  two 
circuits  have  the  same  topology  but  not  necessarily  the  same  circuit 
elements . 
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or 


Vit(s)v;(s)  =^ib(s)v1;(s) 

t  b 


^r(s)vt(s)  -2,vb(8)Ib(8) 


(16) 


(17) 


t  b 

Subtracting  (17)  from  (16)  yields: 


L 

t 


^It(s)V^(s)  -  I«(s)Vt(s)^  =  VLVb(s)Ib(s)  "  Vb(s)Ib(s)]     •   (18) 


Unless  otherwise  stated,  all  calculation  in  the  remainder  of  this 
paper  will  be  in  the  Laplace  domain  (s  -  domain). 
Consider  circuits  as  shown  in  Figure  3. 


For  circuit  (a) 

V(s)  I(s) 

Port      E  (s)  I(s) 
g 

Circuit   V,  (s)  L  (s) 

b  b 


For  circuit  (b) 

V(s)  I(s) 

E  (s)  I(s)  +  61(s) 

8 
Vb(s)  +  6Vb(s)     Ib(s)  +  61  (s) 


Using  equation  (18) , 

I(s)E  (s)  -  [I(s)  +  6I(s)]E  (s)  =X{Ib(s)[Vb(s)  +6Vb(s>3  " 

Vb(s)[Ib(s)  +  6Ib(s)]} 
E  (s)I(s)  =I{vb(s)6Ib(s)  -  Ib(8)6Vb(S)} 


V  (s)  =  Z  (s)I  (s)  and  for  small  changes 
6Vb(s)  =  Zb(s)6Ib(s)  +  Ib(s)6Zb(s) 
Multiplying  (20)  by  I  (s)  and  rearranging  will  yield: 

Ib(s)6Vb(s)  -  Vb(s)6lb(s)  =  Ib(s)5Zb(s) 


(20) 


(21) 


From  the  two  states  in  Figure  3,  it  can  be  seen  that: 

E_(s)  =  I(s)Z(s)  (22) 

and 

g   ' 

=  I(s)Z(s)  +  I(s)6Z(s)  +Z(s)6I(s)  +  6I(s)6Z(s)  (23) 
Substituting  equation  (22)  into  (23)  shows  that: 

I(s)6Z(s)  +Z(s)6l(s)  +  6l(s)5Z(s)  =  0  (24) 


E  (s)  =  [l(s)  +  6l(s)][Z(s)  +  6Z(s)] 

o 
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I  (s)6Z(s)  -  -E  (s)6I(s)  -  I(s)6l(s)6Z(s) 

=  6l(s)[E  (s)  +  I(s)6Z(s)]     .  . 
If  I(s)5Z(s)  «  E  (s)    [i.e.  6Z(s)  «  Z(s)],  then: 

I2(s)6Z(s)  -  -E  (s)5I(s) 

g 

Substituting  (26)  and  (21)  into  (19),  it  is  found  that: 
I2(s)6Z(s)  =Y  I2(s)6Z,(s) 

Lj        D  D 

b 

Since  Z(s)  is  a  function  of  Z  (s) ,  6Z(s)  =  /   "  '  /  N  6Z  (s) 

b  !—>  o^,  (s)    b 

b 
Putting  this  result  into  (27)  results  in 

^Xff^)    5VS>    =Il2(s)6Zb(s) 


Vs)f 


y  az(S)   ._  ,  .    v  r  l 
iszTTs)  6zb(s)  =  Lln^-j  6VS 
b     b  b 


V^2 
.I(s)    J 


(25) 


(26) 


(27) 


(28) 


This  is  Cjhn's  theorem  in  the  s  -  domain  for  any  impedance. 

If  s  is  set  equal  to  jio,  the  results  can  be  obtained  in  the  frequency 
domain : 


zb(jw)      Li(jao)  J 


(29) 


The  right  hand  side  of  (29)  may  be  written  as  a  transfer  function  as  a 
magnitude  and  a  phase  angle. 


^  -  TOO).**"* 
I(ju>)     UJJ)e 


(30) 


sz 


I,  (ju))_2 


(31) 


If  Zb(juD)  =  R, 


giJgl-T2(«o)/29(m) 


BR 


(32) 
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If   Zb(juj)    =   juxL, 


and 


and 


!^   =  T2(.)/2^ 


dZ(jaQ 


SL 


juuT    (uj)/28(oj)    =  wT    (oj)  /26  (a>)    +  f/2 


If  Zb(ju>)    =    1/j'JJC, 


5^)  "  l2«/2iM 


2 
H^-  =    (l/.ju))T2(u))/29(u))    =  ^-^  /29(o»)    -  n/2 


(33) 


(34) 


Thus,    the  magnitude    of    the    sensitivity    changes    as    the    product    of   u>T    (a>) 

2 
or    T    (eu)  /uj   and    its    phase    is    advanced    or    retarded   by   tt/2  when    the    element 

of    interest    is    an    inductor    or    a    capacitor    rather    than   a    resistor. 
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IV.   COIIN'S  THEOREM  IN  THE  TIME  DOMAIN 


Equation  (28)  expresses  Cohn's  theorem  in  the  s  -  domain  for  any 


input  current  I(s)  applied  to  an  impedance  Z(s) 


5Z(s) 
BZb(s) 


I,  (s)n2 

.I(s)  J 


This  equation  is  easy  to  use  if  the  input  current  I,  and  the  internal 

element  current  I,  are  known  as  functions  of  s.   The  input  current  wave 
b 

shape  will  be  known  for  most  application  and  its  transform  to  the  s  -  domain 
will  usually  not  be  difficult.   However,  for  any  but  the  most  elementary 
of  circuits,  the  internal  element  current  I  (s)  can  be  found  only  through 
a  series  of  tedious  calculations.   In  contrast,  there  are  several  good 
computer  programs  which  will  calculate  i  (t) .   It  is,  therefore,  desireable 
to  find  a  direct  means  of  determining  the  sensitivity  as  a  function  of 
time  rather  than  by  taking  the  inverse  transform  of  (28).   For  the  class 
of  input  currents  given  by  (39),  as  derived  later,  it  can  be  shown  that: 

2m 


BZisl  = 
3Zu(s)    L     k 
b      k=0 


k-2n  2,  . 
a.  s     I,(s) 
b 


(35) 


where  n  and  m  are  integers  defined  by  (39).   This  leads  to  the  following 
results  in  the  time  domain: 


Bx 


2n 

=  )    a 


k-2n+j 


b   k=0"  St 


k    k-2n+j  [ib(t)  *  V^ 


(36) 


where  the  symbol  -v  is  used  to  indicate  convolution,  and 
j  -  +1  if  x   is  an  inductor  (L) 

=  0  if  x   is  a  resistor  (R) 

=  -1  if v  is  1/capacitor  (i/C) 
A  very  interesting  simplification  occurs  when  i  (t)  =  6(t),  a  unit  impulse 
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so  that  E(s)  =  Z(s) 


,7  ^iK    s  -  E  (s)I(s) 
BZb(s)/Zb(s)     b     b 

dz,/z,     b      b 
b   b 


(37) 


(38) 


Equation  (37)  shows  that  the  sensitivity  is  related  to  the  instantaneous 
power  in  the  frequency  domain. 

In  attempting  to  develop  the  time  domain  equivalent  of  equation  (35), 
it  would  appear  to  a  casual  observer  that  a  simple  application  of  the 
principal  of  convolution  would  apply.   That  is,  X(s)Y(s)  =  £[x(t)  *y  (t)  ] 
and  thus,  X(s)/Y(s)  =  £[x(t)  >vVy(t)].   However,  this  is  not  true  since 
£[y~  (t)]  r   Y   (s)  .   For  example: 
Let  Y(s)  =  s 

y(t)  =  a  unit  doublet 
Now  let  1/Y(s)  =  1/s 

1 /s  transforms  to  a  unit  step  in  the 
time  domain  and  a  unit  step  is 
not  the  reciprocal  of  the  unit 
doublet  . 
It  would  appear  useful  to  try  to  determine  if  there  is  a  direct  way 
to  obtain  the  inverse  transform  of  1/Y(s)  to  get  a  time  function  which  can 
be  convolved  with  x(t).   There  does  not  appear  to  be  any  such  simple 
operation . 

An  alternate  possibility  is  to  limit  the  consideration  to  input 
current  wave  shapes  which  can  be  manipulated  easily  so  as  to  yield  a  direct 


solution.   Consider  an  input  current  of  the  following  form: 

n 


I(s)  = 


fis 


m 


N(s) 
D(s) 


(39) 


/  b  .s 
l— '     i 
j=0J 
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g  =  any  real  number 

n  =  any  integer  +  ,  -  ,  or  0 

b.  =  all  real   numbers  (all  but  one  can  be  zero). 

This  allows  the  input  to  include  impulses,  steps,  ramps,  sinewaves,  decaying 

.  ..    sin  at   sinh   at 

or  increasing  exponentials,  ,  ,  etc. 

a        a 

Considering  the  input  current  specified  in  equation  (39),  equation  (28) 
becomes 


dZb(5)  N2(s) 


(40) 


2m 

-2    -2n  T""  ,       k   2 

8  s     zAs  Ib(s) 

k=0 

2m 
V  (,       -2    k-2n\2 

k=0 

2m 


■  X(vk-2n>b<s>     •  ^ 

k=0 

The  result  is  equation  (35) . 

Consider  the  impedance  of  the  elements.   If  the  element  is  an 
inductor,  then  Z   =  sL  and  equation  (28)  may  be  written  as 

az(s)    .  .i  (s)  2 

mil .  S(S!^L)2  (A2) 

8LU      \I(s)  /  K      ' 

D 

If  the  element  is  a  capacitor,  then  Z.  =  1/sC   and  equation  (28)  becomes 

b        b 


-I,  (s),2 

r  =  i  \ 

b' 


a(i/c.)     s  vi(s)  ;      •  K^J 


The  subscripts  are  changed  because  the  b's  are  now  the  coefficients 
of  the  squared  polynomial. 
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If  the  element  is  a  resistor,  equation  (28)  is 


5Z(s)  _  /VfV  .       .  t      . 

SRb   -  ^  Ks);     •  (44) 


Thus  in  general  the  sensitivity  can  be  written  as 


Bxb       \I(s)  /  ' 


where  j  =  +1  if  x,  is  an  inductor 
b 

=  0  if  x.  is  a  resistor 
b 

=  -1  if  x,  is  a  compliance 
b 

Combining  (45)  and  (41)  yields  the  result 

2m 
&&L   =  sj   Yfa  sk"2nV2(s) 
b         k=0 


2m 


V  (        k-2n+j\T2,  , 


k=0 
It  can  be  shown  that  multiplication  by  s  in  the  s  -  domain  is  equi- 
valent to  differentiation  with  respect  to  t  in  the  time  domain,  that 
division  by  s  in  the  s  -  domain  is  equivalent • to  integration  in  the  time 
domain,  and  that  multiplication  of  F(s)  in  the  s  -  domain  is  equivalent  to 
convolution  in  the  time  domain.    Thus,  if  the  exponent  of  s  in  a  term  of 
equation  (41)  is  greater  than  zero,  the  individual  terms  of  (46)  become 

k-2n 
at 
If  the  exponent  of  s  is  negative,  the  term  of  (46)  becomes 


See,  for  instance,  Holbrook  [Ref.  4]  pp.  66,  68  and  108. 
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The  number  of  integrations  is  equal  to  k-2n.   The  sensitivity, 


/'BZ 
VoZ, 


then  becomes  a  sum  of  integrations  and/or  differentiations  of  the 
convolution  of  the  current  through  the  element  of  interest  with  itself, 
and  equation  (46)  can  be  written  as: 


az(t)       V       a         j    [" .   ,. v       .   ,.." 

c^~"  =  L-\   ^  k-2n+j  L1b(t)  *  Lb(t). 
b     k=0   dt 


This  is  equation  (36)  . 

To  show  the  power  relations  of  equation  (37),  Ohm's  lav?  is  applied 
(see  Figure  4).   Again,  unless  otherwise  stated,  all  computations  are  made 
in  the  s  -  domain. 

BZ(s)   _   b    J 

sz  (s)      Li  (s)J 

b         g 

Efs)  -  I  (s)Z(s) 
g 


BE(s)  =  I  BZ(s)  +  Z(s)BI  (s) 
g  g 

BE(s) 


az(s)  = 


I  (s) 


Thus 


BE(s)/I  (s)    rIK(s)n2 

£ =  ]_b j 

BZ,  (s)     Li  (s)J 
b  g 


5E(s) 
BZb(s) 


Xb<S> 

1(8) 


BE(s) 


Ib(s)  Ib(s)    Eb(s)Ib(s) 


BZ  (s)/Z  (s)   Z.  (s)  I  (s) 
b     b       b      g 


I  (s) 


(48) 


If  the  input  current  is  limited  to  an  impulse  in  the  time  domain,  I  (s) 
will  be  unity  and  equation  (42)  becomes 


BEXs) 

BZu(s)/Zl(s)   *V  ;Vs; 


(49) 
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If  a  circuit  is  excited  with  an  impulse  current  (see  Figure  5) ,  the 
ratio  of  the  voltage   across  any  element  to  the  input  current  may  be 
written  as  a  transfer  function  such  that: 


I  (s) 


=  T(s) 


(50) 


E(s)  =  T(s)Ib(s) 


If  the  circuit  is  now  excited  by  a  current  source  equal  to  the  current 
through  the  element  of  interest,  the  new  element  voltage,  E(s),  will  be 

(51) 

This  is  true  since  the.  transfer  function  T(s)  is  a  function  of  the  network 
and  not  of  the  input  current  or  voltage.   The  original  input  current  was 
an  impulse  so  that  equation  (50)  can  be  written  as  E  (s)  =  T(s) . 


Substituting  this  into  (51)  will  yield 


E(s)  =  E  (s)I  (s) 
b     b 


(52) 


The  input  current  can  be  written  in  the  time  domain  as  i,  (t)  ,  and  E,  (s)I,  (s) 

b  b     b 

can  be  written  as  e  (t)  *  i  (t)  .   Thus,  (52)  becomes  in  the  time  domain 


e(t)  -  eb(t)  *  ib(t) 


(53) 


Equation  (53)  is  the  right  hand  side  of  equation  (38).   Thus,  in  order  to 

measure  the  sensitivity  of  equation  (38),  the  following  procedure  can  be 

used.   Excite  the  original  circuit  by  an  impulse  and  calculate  i  (t) . 

Then  excite  the  circuit  wi  th   i,  (t)  .   The  voltage  across  the  element  is  then 

b 

equal  to  the  sensitivity. 

There  are  several  interesting  variations  of  equation   (28).   For 
instance:   since  Y(s)  =  Z   (s) 

5Y(s)      = 


i      az(s)  i     rVs)f     r  Vs) 

*VS>     '   Z2(s)      dVS>      '    "z2(s)    Ll<S>    J 


-Ks)Z(s) 


M>)~J 


by(s)    _    rVs)f 

SZb(s)  L  E(s)j 


(55) 
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3Z.  (s)/Z.  (s)   ^b1  ;Le(s)  j      Le(s)  j 

b      b 

If  I  =  EY  is  substituted  into  equation  (28) 

I  (s)  2    rE  (s)Y  (s)  2 
5Z(s)   _  j  b       _    b     b    | 

SZb(s)  "  Ll(s)  J   "  LE(s)  Y(s)  J 

Y2(s)5Z(s)    fEb(s)12 

Y2(s)SZh(s)  "  LE<S>  J      ' 
b      b 

- 1  < 

Remembering  that  Z  =  Y   ,  the  partial  derivative  yields  5Z  =  -cVY/Y 

Putting  this  into  (56)  results  in 

syis)_  _  rVs)f 

BYb(s)    L  E(s)_ 


For  another  variation,  Z  =  Y    is  used.   Thus 


Eb(s),2 


az(s)    =        i       sy(s)   =     i       [_V_n 

aYb(s)      Y2(s)   BYb(s)--Y2(s)  LE(s)  J 


(56) 


(57) 


az(s)    ..  rEb(s)f 

BYb(s)    l  I(s)J 


az(s) 

SY  (s)Y  (s) 
b     b 


(s)E2(s)     Eb(s)Ib(s) 


Ybv"'  ~bv'       "b 

2       '       2 
I  (s)  I  (s) 


(58) 


(59) 


If  e 


(t)  is  restricted  to  a  delta  function,  then  equation  (55),  (56),  and 


(57)  may  be  written  as 


SXi§i    =  -i2(s) 

9Zb(s)     \{S) 

>7    /ds^/gN   =  -E,  (s)I  (s) 
8Z  (s)/Z  (s)       b     b 


(60) 
(61) 


SY(S)  =  E2(s) 
BY,  (s)    Vs; 
b 


(62) 


If  i(t)  is  restricted  to  a  delta,  equation  (58)  and  (59)  may  be  written 


as 
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aVo/vo  -  -V>V>  (64> 

If  the  denominator  of  the  right  hand  side  of  these  equations  are 
restricted  to  the  forms  of  equation  (39),  time  domain  functions  similar  to 
that  of  equation  (47)  can  be  obtained.   For  any  sensitivity  calculation 
however,  the  impulse  driving  function  gives  the  simplest  and  most  direct 
result . 
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Figure    1 

'J\vo  Terminal  Network  with   Variable   r, 

k 

Used  for  Development  of  Cohn's  Theorem 


22 


Figure  2 

Five  Branch  Network  Used  to 
Demonstrate  Colin 's  Theorem. 
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Figure  3 


Two  Circuits  with  the  Same  Topology 
But  with  Different  Internal  Impedances 
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+ 


Figure  4 

Two  Terminal  Networks  Used  for  Developing 
the  Power  Relations . 

If  a  constant  input  current  is  required,  a 
current  source  is  connected  to  the  input 
terminals.   If  a  constant  voltage  source  is 
required,  a  voltage  source,  is  connected. 
In  either  case  E  =  IZ . 
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V->  CD 

i 

1 

K.) 


T(s) 


Ib(s)  Eb(s) 


Ib(s) 


9 

Hh 

t 

E 


E  =  T(s)I  (s) 

b 

"  Vs)Vs) 


£«  a 


Figure;  5 

Two  Terminal  Network  Used  for  Developing 
the  Time  Domain  Solution. 
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APPENDIX  A 
Example  Using  a  Parallel  RLC  Circuit 

-4        -7 
For  the  circuit  in  Figure  6;  let  R  =  1.0,  L  -  10   ,  C  -  10   . 

Calculate  the  sensitivity  of  the  input  impedance  with  respect  to  each  of  the 

circuit  components  for  a  delta  input  current.   Also,  determine  the  sensitivity 

of  the  input  voltage  with  respect  to  6Z/Z  for  each  component  with  the  same 

delta  input . 

An  impulse  was  approximated  by  a  triangular  pulse  8.0(10   )  seconds  wide 
and  2.5(10  )  amps  at  its  center.   This  input  current  was  used  as  data  for 
an  ECAP  solution  to  arrive  at  the  currents  through  and  voltages  across  the 
components . 

The  data  obtained  from  the  ECAP  solution  was  used  as  input  data  to 
Subroutine  SEN  which  is  in  the  Computer  Program  section  of  this  paper.   Since 

outputs  were  obtained  from  ECAP  every  50  calculations,  the  time  increment 

-6  -4 

for  SEN  was  5(10   ).   The  final  time  was  5(10   ).   The  delta  was  chosen 

to  give  it  unit  area  so  COEF  is  1.0.   The  Laplace  transform  of  a  delta  is 

a  constant.   Thus,  SN  and  IPR  were  set  at  zero.   The  IELEM  and  ITYPE>  values 

were  determined  by  the  type  of  element  and  type  of  sensitivity  required. 

The  results  are  shown  in  Figures  7  through  12 . 
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Figure  6 
Parallel  RLC  Circuit  Used  for  Sample  Calculations 
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SUBROUTINE  SEN 
1B,E,T,  SENS) 


(H,TF,COlF,SN,  I  PR,KEY1,I£LEM,  I  TYPE,  A-, 


PURPOSE 

TO  COMPUTi 
OR 


THE  SENSITIVITY 


01  THE  IN?ir  IMPEDANCE 

INPUT  ADMITTANCE  OF  AMY  RLC  CIRCUIT  WITH  RESPECT  TO 

I.  THE  CHANGE  IN  ANY  INTERNAL  IMPEDANCE 

Z.     THE  CHANGE  IN  ANY  INTERNAL  ADMITTANCE 

J.  rHE  CHANGi  IN  ANY  IMP£DANi  "   >IVI0ED  BY  ~HA"r 

[MP!  DANCE 

•fc  Th   •  HAN(    IN  AN  ADMITTANCE  DIVIDED  P>Y  THAT 


! x  M 1 

THE  SUBROUTINE  WILL  ALSO  CALCULATE 
THE  [NPUT  VOLTAGE  wRT  DELTA  7  OV  ■  •< 
INTERNAL  I  MPEDfl  VJCE.  TNT  L  :CTR  ICAL 
EXPRES!  ID  AS  A  LAPL'ACE  TRANSFORM. 


THE  SENSITIVITY  OF 

Z,  WHERE  /  IS  ANY 


SOURC  : 


■*  : 


USAGe 

CALL  St-N(HfTFf(  :EF  ,SN,  I  Pi 


,Ki  Yl  ,IEt  EM,  1  TYPE,  A,  !3,E,SENS  I 


TWO  DAT 4  CARDS  AR, 


REQUI 

C  "r    '. 


THE:  F]  KS  r  *6    CHAPACT! 
CARD  THE  SECOND  V? . 

DESCRIPTION  OF  PARAMETERS 


■i  FOR  THE  GRAPH 
THE  r  IK  ST  C  \RD 


TITLE  IF  A 
list  CONTAIN 


THE  T  ITLE  AND  THE  SECOND 


H      -  THE  INCR<=MFNT  OF  TH-  ARGUMENT  VALUrc 

TF     -  THE  «=INAL  TIME  OF  THE  INPUT  VECTOR  (A) 

COEF   -  VALUE  OF  THE  SOURCE  NUMERATOR  COEFFICIENT 

SN     -  POWER  np  the  S  TERM  Or  TH"  SOURCE  NUMERATOR 

1PR    -  HIGHEST  POWER  IN  THE  SOURCE  NUMERATOR 

KEY]   -  CONTROLS  GRAPH  OUTPUT:  0  PRODUCES  NO  GRAPH 

IELEM  -  IDEMTIFI-S  THE  ELFMFNT  OF  INT^R,-?T:  ENTER  A  1 

FOR  AN  INDUCTOR,  A  0  F^R  A  RESIST  OP,  OP  A  -1 

FOP  A  CAPACITOR 
ITYPE  -  TYPE  OF  SENSITIVITY  REOUIRFD.   ENTER  A  0  IF 

THE  CONVOLUTION  IS  TO  BE  BETWEEN  DIFFERENT 

VECTORS.   ENTFR  fi     1  IF  THE  CONVOLUTION  IS  "n 

BETWEEN  IDENTICAL  VECTORS. 
A      -  CURRENT  OR  VOLTAGE  VECTOR  FOR  TH.E  ELEMENT  OF 

INTERFST 
B      -  CHRPENT  OF  VOLTAGE  VECTOR  FOP  THE  ELEMENT  OF 

INTEREST. 
E      -  COEFFICIENTS  Or  THE  SOURCE  DENOMINATOR  IN 

ORDER  OF  INCREASING  POW-R. 
T      -  OUTPUT  VECTOR  Qe  THE  TIML: 
SENS   -  OUTPUT  VECTOR  Oc  THE  SENSITIVITY 

REMARKS 

1.  THE  EQUIDISTANT  VECTOR  OF  THE  ELEMFNT  CURRENT  OR 
VOLTAGE(A)  MUST  BE  FURNISHED  WHEN  THE  S FNS T T I  VI TY 
SOLUTION  REQUIRES  THE  SELF  CONVOLUTION  OF  CURRENT 
OP  VOLTAGE.  THE  EQUIDISTANT  VECTOR  (B)  MUST  BE 
FURNISHED  ALSO  BUT  m^y  BE  7EP0  UNLESS  THE 
SENSTTIVITY  SOLUTION  REQUIRES. THE  CONVOLUTION  OF 
DIFFERENT  VECTORS.  BOTH  VECTORS  MUST  ^F  IN  5E15.7 
FORMAT. 

2.  MAXIMUM  LENGTH  OE  (A)  ANO  (R)  IS  100  POINTS. 

SUBROUTINES  AND  FUNCTION  SUBPROGRAMS  REOUIRFD 
DAHC02  {H,TF1 ,TF2, A,B,P) 
DET3  (H,Y,7,NDT  »',  IER) 
QSF  {H,Y,Z,NDI«) 
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c 
c 

c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 


c 
c 
c 


c 
c 
c 


c 
c 
c 


c 
c 
c 

c 
c 
c 


METHOD 

SQUARES  A  POLYNOMIAL,  FXTPACTS  THE  POWERS  AMD  USCS 
THESE  POWERS  TO  DETERMINE  r"::    NUMBER  Oc  DIFFERENTIA- 
TIONS OR  INTEGRATIONS  TO  BE  PERFORMED. 

REFERENCES 

(1)  JAMFS  G.  HOLBRO^K,  LADL4CF  TRANSFORMS  EO» 

FLECTRONIC  FNGINFERS,  PERGAMQN  PRESS,  1966. 
(?)  A.L.  DAHLVIG. APPLICATIONS  °r-  COHN*S  SENSITIVITY  ■ 

THEOREM  TO  TIME  DOMAIN  RESPONSES,  ^ASTERS  THESIS, 

NAVAL  POSTGRADUATE  SCHOOL,  1970. 

DI ''ENS  I  ON  A  (100)  iR(100)fC(100f  100)  ,  D  {  2  01  ,  2  )  ,  Y  (  201  )  , 
1Z(  20  1  )  ,W<  2  01)  ,SENS<  201  )  ,E(  11  )  ,DEN<  25)  ,T(?01) 
INTEGER   SD(25),SN 
REAL  LAB1/4H    7 
REALMS  ITITL(12) 

I2=TP/H 

13=12+12+1 

SN=SN+SN 

IP1=IPR+1 

L=IPR+IPR+1 

DETERMINE  THE  TYPE  OF  SENSITIVITY  REQUIRED 


81 


IF  ( ITYPE.FO.O)  GO  TO 
DO  84  1=1,12 
B(I )=A( I) 

84  CONTINUE 
81  CONTINUE 

INITIALIZE  DEN  AND  SO 

SO (3  1=0 

DEN( l)=C.OOO 
DO  1).  1=2,  L 
DEN( I ) =0.000 
SD( I )=I-1 

11  CONTINUE 

SQUARE  THE  INPUT  POLYNOMIAL 

DO  12  1  =  1, 1  PI 
DO  13  K=1,IP1 
J=I+K-1 
DD  =  E(  I  )*E(K) 
DEN( J )=DD+DEN( J) 
1^  CONTINUE 

12  CONTINUE 
COEF=COEF*COEF 
DO  4  1=1, L 

IF(DENU).EQ. 0.0000)00  TO  14 
GO  TO  4 

14  SD(I)  =  0 
4  CONTINLE 

WRITE( 6,906) 

WRITEC  6,111)  (  I,DFN( I  )  ,SD( I  )  ,1=1  ,L) 

PREPARE  FOR  CONVOLUTION 

TFZ=TF 

CONVOLVE  INPUT  VECTORS;  INITALIZE  SENS 

CALL  DAHC02  ( H , TF , TF 7 , A , B , D ) 
DO  16  1=1,12 
T(  I  )=D(  1,1) 
Y(  I  )=D(  1,2) 
SENS(  t  )=0. 00000 
16  CONTINUE 
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c 
c 

c 


c 
c 
c 
c 

c 
c 
c 


DO  21  M=1,L 
J=SD( M)-SN 

DETERMINE  ELEMENT  TYPE:"  AMP  SET  VALUE  HF  J 

IF  UELEM)  74,75,76 

74  J  =  J-1 

GO  TO  75 
76  J=J+1 

DETFPMINF  WHETHER  INTEGRATION,  DIFFERENTIATION  OR 
NO  ACTION  IS  REQUIRED 


75  IF 


22 


(J)  22,23,24 

BEGIN  REQUIRED  INTEGRATIONS  AND  EVALUATE  SENS 


J=—l  Jk  J 

CALL  OSF(HtY,W,I2> 

WRI TE( 6,907) 

IF  (J.LT.2)  GO  TO  73 

DO  31  K=2,J 

CALL  0SF(H,W,W,I2) 

WRITE (6,907) 

31  CONTINUE 
GO  TO  73 

73  DO  4  1  1=1,12 

SENSJ I >=SENS(I )+(DEN(M) /COEF)*W< I ) 

41  CONTINUE 
GO  TO  21 

NO  ACTION  IS  REQUIRED,  EVALUATE  SENS 

SENS( I  )=SENS< I  )+<OEN(M) ^COEF)*Y< I ) 

42  CONTINUE 
GO  TO  21 

BEGIN  REQUIRED  DIFFERENTIATIONS  AND  cVALUATc  SENS 

24  CALL  DET3(H,Y,W,I2,IER) 

WRITE  (6,10M  IER 

I  F( J.LT.2)  GO  TO  99 

DO  3  2  K  =  ?,J 

CALL  DET3(H,W,W,I2 ,IEP  ) 

WRITE  (6, 104}  IER 

32  CONTINUE 

99  DO  43  1=1,12 

SEMS(  I  )=SENS< I  )+(DEN(M) /COrF)*W( I) 

43  CONTINUE 
21  CONTINUE 

DETERMINE  ELEMENT  TYPE  FOR  OUTPUT  PURPOSES 

IF  ( IELEM)15,77,78 

15  WRITE (6,909) 
GO  TO  64 

77  WRITE  (6,910) 
GO  TO  64 

78  WRITE  (6, CO*) 

64  IF  (KEY1.EQ.0)  GO  TO  79 

READ  (  5, 911 )( ITTTL ( I  )  ,1=1,12) 

CALL  DP AW( 12, T, SENS ,0,0, LABI , I TITL, 0,0,1 ,1,1,1,6,4,1, 
1LAST) 

79  CONTINUE 

101  FORMAT (5E15. 71 

104  FORMAT* 5X,1  IER  =  • ,15) 

HI    FORM AT ( 3X , 13 , 3X , E 1 5 . 7 ,6X , I  3 ) 

906  FORMAT ( ?X,  'TERM ■ ,5X, 'COEFF ICI-NT • ,7X, «  P DWEP '  ) 

907  FCRMAT( 1H0,5X, 'SPLIMT' ) 

908  FORMAT ( 1H0,5X , »THE    ELEMENT    IS     AN    INDUCTOR') 

909  FORMAT ( lHOt5X, • THE    ELEMENT    IS    A    CAPACITOR1) 
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910  FORMAT (lH0T5Xt "THE  ELEMENT  IS  A  RESISTOR'} 

911  FORMAT  (668) 
RETURN 

END 


SUB 
DIM 
N=T 

M  =  T 

NM= 

DO 

D(  I 

DO 

DO 

C(  I 

K=I 

D(K 

CON 

D(  1 

DO 

AA  = 

0(L 

RET 

END 


F  DAHC02  (H,TF1  ,TF?f  A,  B,D) 
A(100),B(100),C(100,100),D(201,?) 


ROUT  IN 

ENSIGN 

F  1  /  H 

F2/H 

N  +  M  + 1 

1  I=1,NM 

t2)=0, 

5  I=1,N 

<♦  J  =  1,M 

,J)=A( I)*B( J) 

+  J 

,2) 

TIN 
-1  ) 


=D(K, 2)+C( I, J)*H 
UE 
=  0. 
=  2,NM 


7  L 
L-l 

, 1 ) =  AA*H 
URN 


C 
C 

c 
c 

c 
c 

c 
c 


c 
c 


c 
c 

c 


c 

c 


c 
c 


SUBROUTINE  DFT3 ( H , Y , I , NDI M, IFF) 


DIMENSION  Y(  1)  f  HI) 

TEST  OF  DIMENSION 
IF(NDIM-3)4, 1,1 

TEST  OF  STEPSI7E 

1  IF(H)2v5f2 

PREPARE  DIFFERENTIATION  LOOP 

2  HH=.5/H 
YY=Y(NDIM-2) 
B=Y( 2) +  Y( 2) 

B=HH*( B+B-Y( 3)-Y(  1.  )-Y(  1  )-Y(  1  )  ) 

START    DIFFERENTIAL  CN    LOOP 
DO    3    I=3tNDIM 
A=B 
R=HH*( Y( I )-Y( 1-2) ) 

3  Z( I-2)=A 

END  OF  DIFFERENTIATION  LOOP 

NORMAL  EXIT 
IER  =  0 

A=Y( NDIM-1 ) +Y( NDIM-1) 

Z(NDIM)=HH*(  Y(  NDI  M)+-Y(NDI  W)+Y(  NCI M ) - A-A+YY ) 
Z(NDIM-1)=B 
RETURN 

ERROR  EXIT  IN  CASE  MDIM  IS  LESS  THAN  3 
A  IER=-1 

RETURN 

ERROR    EXIT    IN    CASE    OF    ZERO    STEPS  1 7C 
5     IER=1 
RETURN 
END 
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c 
c 
c 


c 
c 


c 
c 

c 

c 


SUBROUTINE  OSF( H , Y , 7. , NDI M) 


DIMENSION  Y(  1)  ,7.(1) 

HT=.3333333*H 
[F(NDIM- 5)7,8,1 

NDIM  IS  GREATER  THAN  5.  PREPARATIONS  OF  INTEGRATION 
LOOP 

1  SUM1=Y(2)+Y(?) 
SUM1  =  SUNT1  +  SUM1 
SUM1=HT*(  Y(  1  )*S'.IMH-Y(3)  ) 
AUX1=Y(4)+Y(4) 
AUX1=ALX1+AUX1 

AUX1  =  C.IJV1+HT*(  Y(3)  +AUX1  +  Y(  =  )  ) 

AU X2=HT*< Y( ] )+3. 87  5 *(Y(?)+Y(C}) +2.625 *(Y(3)+Y(4))+ 
1Y(6) } 
SUM2=Y( 5) +  Y( 5) 
SUM2=SUM2+SUM2 

SUM2  =  AUX2-HT*(  Y(4)  *$I|M2+Y(M  ) 
7.  (  1 )  =  0  . 
AUX=Y( 3)+Y(3) 
AUX=AUX*-AUX 

7(  2)  =SUM2-HT*<  Y(  2)  +  AUX+YU  )  ) 
7( 3)=SUM1 
7(4)=SUM2 
IF(NDIM-6)5,5,2 

INTEGRATION    LOOP 

2  DO    4    I =7, NDIM, 2 
SUM1=AUX1 
SUM2=AUX2 

AUX1=Y(  I-l)+Y(  1-1 ) 
AUX1  =  AUXH  UJX1 

AUXl  =  SUMH-HT*(Y(I-2)+AUXl  +  Y(in 

7(  I-2)=SUM1 

IF(  I- NDIM)  3,^,6 

3  AUX2=Y( I )+Y( I) 
AUX2=A  LX2  +  aiJX? 

AUX2=SLM2+HT- ( Y ( 1-1  )+AUX2  +  Y< 1+1 ) ) 

4  Z(I-1)=SUM2 

5  Z(NDIM-l )=AUX1 
Z(NOIM)=AUX2 
RETURN 

6  7(NDIM-1)=SUV2 
Z(NDIM)=AUX1 
RETURN 

END  OF  INTEGRATION  LOOP 

7  IF(NDIM-3)12,11,8 

NDIM    I S    EQUAL    TO    4    OR    5 
P    SUM2=1.125*HT*(  Y(1)+Y{2)+Y(2)+Y<2)+Y(3>*Y(?H-Y{3)  + 
1Y(4) } 

SUM1=Y( 2)+Y(2> 

SUM1=SUM1+SUM1 

SUM1=HT*< Y( 1 )+SU*U  +  Y(  3 )  ) 

7(  1  )=0. 

AUX!=Y( 3)+Y(  3) 

AUXl=AtXl+AUXl 

7(  2)=SUM2-HT*< Y( 2)+AUXl+Y(4)  ) 

IF( NDI M-5) 10,9,9 
9    AUX1  =  Y(4H-Y(4) 

AUX1=ALX1 +  AUXI 

Z(5)=SUM1+HT*( Y(3)+AUX1+Y(5)) 
10    Z<3)=SUM1 

Z(4)=SUM2 
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c 

c 


RETURN 

NDIM     IS    EQUAL    TH    3 
li    SUM1=HT*( 1.25#Y( 1) +Y(2)+Y( 2)-.2  5*Y (3) ) 

SUM2=Y(  2)+Y(2) 

SUM2=SUM2+SUM2 

Z(  3)-HT*(  Y(  1)+SUV2+Y(3)  ) 

Z(  1)=0. 

Z(2)=SUM1 
12    RETURN 

END 
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